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NOMENCLATURE 


C molecular velocity vector 

d distance between the plates 

f velocity distribution function 
fi , f 2 components of two-stream Maxwellian 

F external body force 

g, G functions defined by equation (45) 

l n (p) integral defined by equation (3) 

f, j, k unit vectors in the x, y and z directions, respectively 
k thermal conductivity 

L plate temperature ratio parameter, equation (43) 

m mass of molecule 

n molecular number density and parameter in equation (3) 
n -\ , n 2 number density functions 

n unit normal vector to an area element 
p pressure 

q heat flux vector 
q 2 heat flux in y direction 

Q, \| / molecular property 
R gas constant 
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T i , T 2 temperature functions 

T y , T L upper and lower plate temperatures 

u, v, w components of p in the x, y and z directions, respectively 

V macroscopic velocity vector 

(Xj, aj constants, equations(51--53) 

P constant in equation (3), and parameter defined by equation (54) 
7 the ratio of specific heats 

« azimuthal angle in binary collision 
n dynamic viscosity 

p mass density 
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SUMMARY 


Plane compressible Poiseuille flow with uniform pressure (Couette flow with stationary 
boundaries) is revisited where the Lees two-stream method with the Enskog equation of 
change is applied. Single particle velocity distribution functions are chosen, which pre- 
serve the essential physical features of this flow with arbitrary but uniform plate tempera- 
tures and gas pressure. Lower moments are shown to lead to expressions for the 
parameter functions, molecular number densities, and temperatures which are entirely in 
agreement with those obtained in the analysis of Lees for compressible plane Couette 
flow in the limit of low Mach number and vanishing mean gas velocity. Important simplifi- 
cations result, which are helpful in gaining insight into the power of kinetic theory in fluid 
mechanics. The temperature distribution, heat flux, as well as density, are completely de- 
termined for the whole range of Knudson numbers from free molecular flow to the contin- 
uum regime, when the pressure level is specified. 



INTRODUCTION 


The Kinetic theory of gases is one of the most intriguing inventions in the theory of matter, 
with its foundations in statistical physics and applications in fundamental fluid physics. 
The main success of the kinetic theory is explaining the transport properties of gases set 
forth by its founders, Clasius, Maxwell (ref. 1 ) and Boltzmann (refs. 2 and 3), and later by 
the more exact works of Chapman and Enskog (ref.4). There are as yet unanswered in- 
triguing questions, -for example, why the Navier-Stokes equations describe fluid mo- 
tions so satisfactorily. They can be obtained from the Boltzmann equations by the 
Chapman-Enskog theory by successive approximation with m as the expansion param- 
eter and a acts on the five macroscopic variables, density p, velocity V, and temperature 
T, where x is the mean free path. Second order hydrodynamical equations derived by Bur- 
nett, however, lead to unrealistic physical phenomenon, -for example, negative absorp- 
tion coefficient in very high frequency sound waves. Turbulence, of course, remains an 
important unsolved problem in fluid mechanics. Attempts by several researchers have 
been made to develop a kinetic theory for turbulence. Despite their modest success, de- 
scription of turbulent flow must clearly be statistical and this makes kinetic theory an at- 
tractive guiding light. 

From the foregoing it may be clear that much insight into the physics of fluids, not only 
turbulence but, for example, the structure of shock waves can be obtained by carefully 
studying the works of masters and understanding their logic in deriving fundamental con- 
cepts and equations and expanding upon them. Even seemingly simple flow geometries 
may serve this purpose enormously. Obtaining exact solutions of Boltzmann equations is 
undoubtedly a formidable task. Lees (ref.5) developed a method utilizing a distribution 
function that is two-sided in character and contains parameter functions of space and time 
variables. The method is formulated in terms of Maxwellian functions which is a natural 
extension and generalization of Mott-Smith's (ref.6) approach for normal shock waves. 
Sufficient numbers of moments are taken to ensure a complete set of first order partial 
differential equations to solve for the functions. As pointed out by Lees (ref.5), this method 
amounts to satisfying the Maxwell-Boltzmann equation in a certain average sense rather 
than point by point. The distribution function is then a suitable weighting function and not 
an exact solution; therefore, the rigidity inherent in the polynomial of Chapman-Enskog 
or Grad’s moment method (ref.7) is avoided. Kinetic theory is often obscured by mathe- 
matical formalism. In this paper we start with a brief discussion of the fundamental con- 
cepts and the two-stream method and proceed with the solution of the present problem in 
a most straightforward way. 
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FUNDAMENTAL CONCEPTS 


Maxwell and Enskog Equations of Change 

The Boltzmann equation is fundamental to the kinetic theory. It represents conservation 
of molecules in the (c, r) phase space and accounts for all the mechanisms which can 
change the number of particles in (dc, dr). Here, the molecular velocity is denoted by c 

and r is the position vector. The RHS of the Boltzmann equation pertains to gas phase 
binary collisions, where no correlation between particle velocities prior to collisions exists. 
This is generally referred to as independent a priori probabilities or molecular chaos. A 
spherically symmetric intermolecular force field which is much stronger than all external 
forces is assumed . The duration of collisions is also much smaller than the time between 
collisions. In an integral form the Boltzmann equation leads to the Maxwell equation of 
change 


|jQfdc + ||;- Jofcdc- Jff ~dc = JJJJ Q “ ff 1 ) gbdbdedcdc' (1) 


c'ceb 


ri * 

where, ^ represents the partial derivative with respect to time. Here F denotes the exter- 
nal force field and e is the azimuthal angle in the collision. The distribution functions f and 
fi represent the binary (probe and colliding) particles before the collision and the primed 
functions after the collision. The speed of approach of binary particles is denoted by g, 
and b is an impact parameter (a hypothetical distance of closest approach if there was no 

interacting force.) The molecular property Q = Q( c) only. In its more general form, when 
Q = Q(C, r, t), where t is the time, the Enskog equation of change is obtained. 

^jQfdc-h^;- Jofcdc-Jf^ +c-^ + F-^ jdc = Jq (f’f’., - f^) gbdbdedcdc’ (2) 
c 9 9 


Averaging Over the Velocity Space 

Since the transport coefficients p, p, k (density, viscosity, and thermal conductivity, re- 
spectively) are given in terms of the averages over the velocity space involving the distri- 
bution function, the conservation laws can be obtained when Qs are the collisional 

invariants, mass m, linear momentum me, and translational kinetic energy of particles 
1/2 me 2 . The RHS of equation (2) is then identically zero. 

In evaluating the three-fold integrals over the velocity space, since the three components 
of the particle velocity are independent, the integrals can be separated into single-fold in- 
tegrals. The single-fold integrals can then be evaluated individually. We will present sum- 
maries of the integral results, employing the following definition 
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oo 


( 3 ) 


For n = 0, 1: 



For n>1 , a recursion formula may be written 

■n<»> = C^ln-sW 


(4) 

( 5 ) 


(6) 


Lees Two-Stream Method 

The integral equations of transfer, equations (1) and (2) possess an important property. 
They permit a large degree of flexibility in the choice of the distribution function. The dis- 
tribution function f can be expressed in terms of arbitrary functions of space and time 
that preserve the essential physical features of a problem and are, therefore, tailored to 
the particular case. There is no sequence of logic which enables one to find such func- 
tions; the tools seem to be accumulated experience and intuition. As mentioned earlier, 
Lees (ref. 5) developed a method utilizing a distribution function which is two-sided in 
character and contains parameter functions of space and time variables. The method is 
formulated in terms of Maxwellian functions, which is a natural extension and generaliza- 
tion of Mott-Smith's (ref. 7) approach for normal shock waves. Sufficient numbers of mo- 
ments are taken to obtain a complete set of first order differential equations to solve for 
the functions. As pointed out earlier, this procedure amounts to satisfying the Maxwell-- 
Boltzmann integral equation and, therefore, the Enskog equation of change in an average 
sense rather than point by point. The search for higher order macroscopic equations in 
terms of the mean quantities is not persued, but surface boundary conditions can be in- 
corporated. Lees (ref. 5) applied this method to steady plane Couette flow for the limiting 
case of low Mach number and low temperature difference. Reasonable results were ob- 
tained for macroscopic quantities over the whole range of densities from free molecular 
flow to the Navier-Stokes. Later, Lees and Liu (ref. 8) expanded this work to compressible 
plane Couette flow with arbitrary Mach number and uniform plate temperature. 

We will start with a two-sided velocity distribution function, which appears sensible for the 
present boundary conditions and carry out the integrations over the corresponding veloc- 
ity space. The results are then shown to be entirely consistent with those obtained in the 
more general analysis of Lees and Liu (ref. 8) for compressible plane Couette flow in the 
limit of low Mach number. Important simplifications result which are helpful in gaining in- 
sight into the power of kinetic theory in fluid mechanics. The temperature distribution, heat 
flux, as well as density, are therefore completely determined for the whole range of Knud- 
son numbers from free molecular flow to the continuum regime when the pressure level 
is specified. 
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PLANE POISEUILLE AND COUETTE FLOWS 

Plain Poiseuille flow and plain Couette flow must be represented by different velocity dis- 
tribution functions in the kinetic theory. They are, however, physically identical when, in 
each flow respectively, uniform gas pressure and stationary boundaries are approached. 
We consider two stationary infinite parallel plates designated by the subscripts U (upper) 
and L (lower) which are held a distance d apart and are at temperatures Ty and T|_ , 
respectively. We will assume that the gas between the plates is monatomic and is at a 

uniform pressure p . Despite the non-zero molecular velocity C, no shear stresses exist, 
and determination of the heat flux and temperature distribution is of primary interest. 

We choose a two-sided Maxwellian particle velocity distribution function f and assume 
that the plates scatter molecules diffusely with complete thermal accommodation. 


n-. -c 2 /2RT 1 

372 e 

(27CRT,) 372 


, v < 0, all u and w 


(7) 


f 2 


n 2 -c 2 /2rt ? 



(2itRT,) 3/2 


, v > 0, all u and w 


(8) 


where n^ n 2 , Tj , and T 2 are parameter functions of y, and R is the gas constant. The 
boundary conditions for completely diffuse re-emission are 


Ti=Tu 

at y= +1/2 d 

(9) 

t 2 =t l 

aty= -1/2 d 

(10) 


The coordinate system denotes the components of velocity by u, v, and w in the x, y, and 

AAA 

z (i, ], k unit vectors) directions, respectively. The plates are in the x, z plane. 


The Heat Flux Vector 

The net flux of any molecular property X|/ (C) across an area element with unit normal vec- 
tor n may be represented by 

J(n-c)\|/(c)f(c,r,t)dc (11) 

c 

The heat flux vector between the plates is directed perpendicular to the plates, i.e., 
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(12) 


q = q 2 J 

Here f (£) is the molecular kinetic energy 1/2 me 2 . Therefore, 

q 2 = (13) 

c 

Introducing the two-stream distribution function f and integrating over the velocity space, 


q 2 = 


mn 1 
2 (2jiRT 1 ) 


oo Q oo 

372 J J Jv(U 2 + V 2 + W 2 ) 

— OO — OO — OO 


-(u 2 + v 2 + w 2 ) /2RT 

e 


i 


dudvdw 


+ 


oo oo oo 



- (u 2 + v 2 + w 2 ) /2RT ? 

e 


dudvdw 


(14) 


— mn 1 mn 2 

q2= + , + 2 (2 , RT2 )3/2 2 < 15 > 


Here, subscripts 1 and 2 to the square brackets refer to the integrals evaluated by equa- 
tions (3)— (5) containing the respective parameter functions T-j and 7% . Therefore, 



We note that q 2 is independent of y and therefore 

.\-n 1 t| /2 + n 2 T| /2 = const ( 17 ) 

The Temperature and Density Distributions 

We now proceed with the definition of temperature as the mean molecular kinetic energy, 
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|mRT = < > 

2 2 


( 18 ) 


where 


mC . a ■ i i w r _l 

<^r-> = — I — = — fdC 


1 f mc2 
nJ 2 

c 


(19) 


T = 3nR J C V° + J 


c 2 f 2 dc 


(c, v < 0) 


(c, v > 0) 


( 20 ) 


T= 


(2tiRT 1 ) 


3/2 


Ml 


(2jcRT 2 ) 


3/2 


_ 4I 2 I o]< 


(21) 


T = 


Ti n i + T 2 n 2 
2n 


(22) 


n = Jfdc 


(23) 


n = 


n. 


(2ttRT 1 ) (c, v < 0) 


I 


-cV2RT 


dc + 


n. 


(2itRT 2 ) 3/2 (CjV J >0) 


J 


-CV2RT, 


dC (24) 


4n. 


n= 


(2tiRT 1 ) 


3/2 


nsi 1 + 


4n, 


(27iRT 2 ) 


3/2 W 2 


(25) 


.*.n = 


= 7> ( n i +n 2 ) 


(26) 


The temperature T is, therefore, 


.\T = 


Ti n i +T 2 n 2 
n., +n 2 


(27) 


9 



The Moments 


Moments may now be taken to establish additional equations for determination of the 
functions n-| , n 2 ,T-| , and T 2 . Since there is no net flux of mass in the y direction at any ’y’, 


J vf 1 dc+ J vf 2 de = 0 

(c,v<0) (c,v>0) 


-4n 


(2tt:RT 1 ) 


7372 [l 0 ! ll 1 + 


4n 2 

(2tcRT 2 ) 


3/2 


[loM 


= o 


(28) 

(29) 


~ "2/^ (30) 

This is the first moment, or a statement of the conservation of mass. The second and forth 
moments related to the x and z components of velocity are not needed. 

The Enskog equation of change (eq.(2) ) for this problem, when no external forces are 
present is reduced to 


i|Qvfdc-Jfv^ = J Q (31) 

C C 

where only the derivatives with respect to y are retained, and Jq denotes the collision in- 
tegral. 

For the third moment (Q = mv) , the conservation of momentum in the y direction, the col- 

HO 

lision integral vanishes. Furthermore, ^ = o , since v is independent of y. Therefore, 

v 2 f 1 dc+ J v^gdc = const (32) 

(c,v<0) (c,v>0) 


4n 


(27iRT 1 ) 


3/2 ^0*2^ -j + 


4n 2 

(27t:RT 2 ) 


3/2 


iw, = const 


(33) 


.%n 1 T 1 + n 2 T 2 = const (34) 

For the fifth moment, Q = mc 2 /2 , the RHS of the equation of change is zero, since it is 
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also a collisional invariant, and S!P is also zero. 

dy 


= 0 (35, 

c 

We have already evaluated this integral in determining the heat flux and shown that equa- 
tion (1 7) is satisfied. 

The sixth moment, Q = mvc 2 , is not a collisional invariant and evaluation of the collision 

do 

integral is required, although again vanishes. We have 

dy 



J 

mvc 


2 


(36) 


where J 2 denotes the collision integral, 
mvc 

For molecules obeying an inverse fifth power intermolecular force law, the collision inte- 
grals in a full-range moment method are completely independent of the form of the distri- 
bution function. This advantage is preserved in the half-range moment method, at least 
for certain moments. Following Maxwell (ref. 1 ) and Jeans (ref. 9), Lees (ref. 5) considered 
collision dynamics between particles by means of vector geometry and showed that the 
first nine moments beyond the collisional invariants are, in fact, completely independent 
of the form of the distribution function. Lees then evaluated and corrected once and for 

all the numerical errors in the expressions given by Jeans . The result for Q = mvc 2 is 
simply 


mvc 


2 


pq 2 

6jut 


(37) 


where p is the gas pressure and |i the classical coefficient of viscosity. 

As pointed out earlier, q 2 is a constant, therefore J 2 > s a constant. Again, the integral 

mvc 

in equation (36) is evaluated using the two-stream distribution function f| and f 2 . 

J (2v 2 u 2 + v 4 )f 1 dc-i-^- J (2v 2 u 2 + v 4 ) f 2 dc = const (38) 

^ (C, V < 0) (C, V > 0) 
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d 

dy 


1 r..2. , .,2, , n 2 rn.2. _i_ ai 2 , 


L (2jcRTj) 


3^ [8I 2 I 0 + 4I^I 4 ] i + 


(2tiRT 2 ) 


^[sisio+^g 


= const 


''c^[n 1 Ti + n 2 T 2 ] = const 


(39) 

(40) 


Agreement With The Results of Lees and Liu 


Lees and Liu (ref. 8) expanded the analysis of Lees (ref. 5) and applied the two-stream 
method to plane compressible Couette flow with arbitrary plate temperatures. Mach and 
Reynolds numbers were defined in terms of the characteristic quantities n|_ and T[_ (den- 
sity and temperature at the lower plate surface respectively), U the relative plate velocity, 
and the separation distance d. 



(41) 


Re 


mn L Ud 

^L 


(42) 


Classical fluid viscosity coefficient fi L evaluated at the lower plate temperature, the ratio 
of specific heats y, and a plate temperature ratio parameter L were also employed. 


L - a/T ij/T l (43) 

It was shown that when the Mach number is small there is a split in the system of equa- 
tions containing the thermodynamic variables n 1( n 2 ,T 1 , and T 2 from the dynamical vari- 
ables u-j and u 2 which were Included in their analysis for moving plates. In the present 
analysis the Mach and Reynolds numbers are vanishingly small. Their ratio Re/M, how- 
ever, is finite and is proportional to the ratio of the mean free path X to the characteristic 
length d. It characterizes the density level of the gas. We note that there are close simi- 
larities between the results of Lees and Liu (ref. 8) and the present analysis when M « 1 . 

As mentioned earlier, Plain Poiseuille flow and plain Couette flow must be represented by 
different velocity distribution functions in kinetic theory. They are, however, physically 
identical when in each flow respectively, uniform gas pressure and stationary boundaries 
are approached. The present analysis must, therefore, be in agreement with the results 
of Lees and Liu (ref. 8) for M « 1. This, in fact, is the case. We note that by substituting 
equation (30) into equation (1 7), one obtains the following equation. 
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(44) 


n l/M T 2- T l) = const 

Inspection of reference 8 reveals, that equations (22a), (22c), and (22d) of reference (8) 
are identical to equations (30), (34), and (44) of the present work, respectively, except for 
the dimensionality difference. A new function g(y) may now be conveniently defined as 

g (y) = i/(n t /f7) = i/(n 2 /Q ( 45 ) 

As shown in reference 8, the present equations (30), (34), (44), and (45) can be easily 
manipulated to result in relations for Ti,T 2 , n-|, and n 2 in terms of g(y) and constants aj . 
For the present problem a-| =0. 


T i (y) = ^ ( 329 - 33 ) 2 

(46) 

T 2 (y) = i(a 2 g + a 3 ) 2 

(47) 

0l <y) 9 (a 2 g - a 3 ) 

(48) 

n2(y) ” g(a 2 g + a 3 ) 

(49) 


Once g(y) is determined, Ti,T 2 , n-j , and n 2 are also completely determined. In order to ob- 
tain g(y) explicitly, the expressions for Ti,T 2 , n-j , and n 2 are substituted in equation(40) 
and the integration is carried out, which results in an additional constant a* The boundary 
conditions (eqs. (9) and (1 0) ), as well as equations (46)-(49), are fully applicable for de- 
termination of the constants. The details of the algebra are rather lengthy and will not be 
presented here. We will, for completeness, however, include the expressions for the a; 
and g(x) as presented in reference 8. For consistency with reference 8, aj can be replaced 
by 0 Cj, and g(x) by G(x). This implies that the quantities n 1f n 2 ,T -j , and T 2 will be accord- 
ingly dimensionless (normalized by n|_ and T|_, respectively). 


G(y) = 


( / \2 , , n/ x / 32 [ 2 a 3 „ Re .v 

(a 3 /a 2 ) + (2/a 2 ) X ¥> ) 


1/2 


-> 1/2 


(50) 


1 /2 

[(1 +P) (L 4 + 2L 3 + p) L 2 ] - (L 3 + 2L 2 — p — 1 ) 


a 2 = 


1 + P - L? 


(51) 



(52) 


a„ = 2 — a r 


a 4 = 1 - 2(1 + p) X (a 3 /a 2 ) + (a g /.a 2 ) 2 (53) 

where 

P = (8/15) Jz/ny Re/M (54) 


The results of both analysis are entirely consistent with each Other. As shown by Lees 
and Liu (ref. 8), the heat flux and temperature distribution are, therefore, fully determined 
for the whole range of Knudson numbers from free molecular flow to the continuum re- 
gime, when the pressure level is specified. Furthermore, when Ui = u 2 = 0 in Lees and Liu 
(ref. 8), the distribution functions lead to identical expressions. The present analysis, 
therefore, may be viewed as an alternative to the more general analysis of Lees and Liu 
(ref. 8) in which identical results are arrived at independently . 
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CONCLUDING REMARKS 

Plane compressible Poiseuille flow with uniform pressure (which is, in fact, Couette flow 
with stationary boundaries) was chosen to discuss the mathematical and physical fea- 
tures of an approximate technique for the solution of the Boltzmann equation. Lees two- 
stream method, a clever invention in kinetic theory, shows that in a rigorous field such as 
statistical physics, fundamental and physically realistic treatments can be made with rel- 
ative ease. The Enskog equation of change, which is more general than the Maxwell in- 
tegral equation, was shown to lead to results that completely describe the heat transfer 
process. A simple form of the two-stream single particle velocity distribution function 
which preserves the essential physical features of Poiseuille flow with arbitrary but uni- 
form plate temperatures and gas pressure was chosen. Important simplifications result 
which are helpful in gaining insight into the power of kinetic theory in fluid mechanics. Low- 
er moments are shown to lead to equations for four parameter functions, molecular num- 
ber densities, and temperatures employed in the distribution function, which are entirely 
consistent with those obtained by Lees and Liu (ref. 8) in the limit of low Mach number and 
vanishing mean gas velocity. The temperature distribution, heat flux, as well as density, 
are therefore completely determined for the whole range of Knudson numbers from free 
molecular flow to the continuum regime, when the pressure level is specified. 
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